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Abstract 

We studied the pairing instabilities in KxFe2-y^^2 using a two stage functional renor- 
malization group (FRG) method. Our results suggest the leading and subleading pairing 
symmetries are nodeless d^2_y2 and nodal extended s respectively. In addition, despite 
having no Fermi surfaces we find the buried hole bands make important contributions to 
the final effective interaction. From the bandstructure, spin susceptibility and the FRG re- 
sults we conclude that the low energy effective interaction in K^.Fe2-ySe2 is well described 
by a Ji — J2 model with dominant nearest-neighbor antiferromagnetic interaction Ji (at 
least as far as the superconducting pairing is concerned). In the end we briefly mention 
several obvious experiments to test whether the pairing symmetry is indeed dj.2_y2. 

1 Introduction 

Very recently a new wave of excitements occurred in the field of iron-based superconductors. 
This is stirred up by the discovery of Ki.Fe2_j/Se2[lJ. These compounds are isostructural to the 
"122"-family of iron pnictides, e.g. BaFe2As2, with the highest transition temperature ~ 33K 
among iron chalcogenides. The reason these compounds attracted considerable interests is not 
their "high Tc". Rather it is because they have a very different electronic structure from all other 
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iron-based superconductors. In particular, the hole pockets near the Brillouin zone center in all 
other iron-based superconductors are found to be absent in the highest Tc K^Fe2-j,Se2[l2l[3l. The 
result of ref. [j2i suggests a very uniform lOmeV) superconducting gap around the electron 
pockets. 

The absence of hole pockets is expected from valence count: KFe2Se2 should have 0.5 
doped electron per Fe relative to the parent iron pnictides. In view of the wide spread belief that 
the scattering between the hole and electron pockets is crucial for the high pairing scale in the 
iron-based superconductors [ffl |5l[a|7l[i|9l[ia[IIl[I2[Il[Il[l5l[l6l[l7l[ia, itis surprising 
that a material without hole pocket can support such a high Tc. It is fair to say the relatively 
high pairing scale in the absence of hole pockets calls for a re-evaluation of the spin fluctuation 
pairing mechanism. This is so because in the other extreme KFe2As2, which has 0.5 doped 
holes per Fe and only hole pockets, is a very low (~ 3-R") [[T9l l20l superconductor with 
experimental evidences of gap nodes [[2T|. 

In this paper we apply the functional renormalization group (FRG) method[l6l|9l[T4l to study 
the pairing instability of K2;Fe2Se2. When the electron-electron interaction is weak compared 
with the bandwidth, this method is unbiased. It sums all virtual one-loop scattering processes 
including particle-hole, particle-particle and vertex corrections. As other iron-based supercon- 
ductors, the strength of electron-electron interaction in K^Fe2-ySe2 is uncertain. Ideally, we 
should combine FRG with the variational Monte-Carlo calculation [|22l . which is underway. 

2 Model 

According to an earlier DFT result[|23l, KFe2Se2 has cylindrical electron pockets around (vr, 0, k^) 
and (0, vr, k^) as well as a 3D electron pocket centered around (0, 0, vr). In our study we uses a 
two dimensional(2D) five-band tight binding model to describe the A;^ = plane of the band- 
structure. It is obtained via the maximally localized Wannier function fit [|24l to our own DFT 
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calculation. In our DFT calculation the plane wave basis method was used ll25l . We adopted 
the generalized gradient approximation (GGA) of Perdew-Burke-Emzerhof formula [|26l for 
the exchange-correlation potentials. The ultrasoft pseudopotentials [|27l were used to model 
the electron-ion interactions. After the full convergence test, the kinetic energy cut-off and 
the charge density cut-off of the plane wave basis were chosen to be 800eV and 6400 eV, re- 
spectively. The Gaussian broadening technique was used and a mesh of 18 x 18 x 9 k-points 
were sampled for the Brillouin-zone integration. In the calculations, the experimental tetrag- 
onal lattice parameters were adopted, and the internal atomic coordinates within the cell were 
determined by the energy minimization. 

Our bandstructure and the associated Fermi surfaces are shown in fig. [T] Because experi- 
mentally there is uncertainty about the Fe and K content, we have adjusted the chemical poten- 
tial so that the distance from the top of hole bands at T to the Fermi energy is about O.leV. It 
mimics the ARPES finding of ref. Q. The parameters of the tight-binding model involves the 
nearest- and second- neighbor hopping among the Fe 3d orbitals: 

i^band = 5^ XI 5Z ^r/cLS/3s- 

ij aP s 

Here j labels the Fe sites, a, (3 label the five different Fe orbitals {a = 1, ..5 denotes Z'^, XZ, 
YZ, X"^ — F^, XY , respectively), s labels spin. The label X, Y refer to the diagonal directions, 
while x, y refer to the nearest neighbor Fe-Fe directions. 
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Figure 1 : (Color online) The band structure and Fermi surfaces. Note that all the Fermi surfaces 
are electron like. 
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In unit of eV the tight-binding parameter are given by 



/ 7.13365 
7.38168 



,a/3 








7.31272 











a/3 



' 



t: 



i,i+y 



7.3816 



\ 7.06574 / 

/ 0.00376 0.07892 0.07892 -0.25696 \ 

0.07892 -0.1406 -0.05034 0.18873 0.20995 

0.07892 -0.05034 -0.1406 -0.18873 0.20995 

0.18873 -0.18873 -0.10725 

\ -0.25696 0.20995 0.20995 -0.35102 / 

/ -0.00559 -0.13684 0.11502 \ 

0.07104 0.05964 

0.13684 0.23042 0.00961 

0.11502 -0.00961 0.08569 

V -0.05964 -0.11843 / 

In terms of these parameters the 5x5 Bloch Hamiltonian is given by 
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Figure 2: (Color online) The real(left) and imaginary(right) part of the bare spin susceptibility. 

are the matrices representing the combined 90° rotation and z reflection and the combined z 
reflection and a gauge transformation on one Fe sublattice, respectively. 

The bare magnetic susceptibility computed using the above bandstructure is shown in fig.|2l 
A broad maximum around Q = (vr, tt) is found in Re(x). If we make an analogy with the 
pnictide superconductors, where Q = (tt, 0) or (0, vr) and the predicted pairing form factor is 
s±, it is natural to guess the pairing form factor here to be (i2,2_j^2[|3l |28l. Note that this gives 
rise to fully gapped superconducting state with the order parameter changing sign between the 
electron Fermi surfaces. In the following we check whether this is true via a FRG calculation. 

As in ref. [|9l we model the electronic correlations by the Hubbard and Hunds types of local 
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interactions. The Hamiltonian is given as follows: 

H =i/band + Ui ^ Uia^^riial + U2 ^ niaUiii 
ia i,a<f5 

+ [ XI 5Z ^^c.<r4f3a' Ciaa' Cip„ (9) 
Q</3 era' 

+ 4at4aiCil3iCi,3^ + h.C. . 

In the rest of the paper we use (f/i, f/2, </h)=(4, 2, 0.7) eV respectively. These value are 
suggested by early dynamical mean-field theory results for the pnictides l|29ll . The sole rea- 
son for using this set of parameter is they represent intermediate couplings, which we believe 
K^Fe2_2;Se2 is likely to be. 

3 FRG Method 

Details of the FRG method are discussed in ref. [I9l[l4l. We made one important modification to 
include the buried hole bands in the FRG calculation. In our previous treatment bands without 
Fermi surface were ignored completely. In the case of Ki:Fe2Se2 (and possibly in heavily hole- 
doped KFe2As2 as well) we believe this is no longer legitimate. The reason is two-fold. First 
consider the three different energy scales in the iron-based superconductors, the largest energy 
scale is the bare bandwidth and bare Coulomb interactions of order eV, the second energy scale 
is the effective spin exchange interactions about 50meV as indicated by neutron studies [l30l . 
the pairing scale(gap) is the smallest about ~ lOmeV. The distance between the top of the 
hole bands and the Fermi level in Ka;Fe2Se2 is smaller than the first but comparable to the 
second energy scale. Therefore for a wide range of cutoff energies virtual excitations from the 
hole band to above the Fermi level should contribute to the renormalization of the effective 
interaction. The second reason is that a preliminary FRG which only include the single band 
that intersects the Fermi energy did not produce any strong pairing tendency, except very weak 
Kohn-Luttinger-type pairing in high angular momentum channel. 
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Physically when the FRG cutoff energy A is much larger than the distance AEq between 
the top of the hole bands and the Fermi energy, the fact that there are no hole pockets should 
not be important. Only when A < AEq the lack of hole pockets will play an important role. 
Consequently, our FRG scheme consists of two stages. In the first stage the FRG scheme is the 
same as the case where the hole pockets are present. Thus we discretize the electron pocket 
band into patches represented by Fermi surfaces points as before. For the hole bands we put 
fictitious "Fermi surfaces", small circles, around (0, 0) and (vr, tt) and discretize the Brillouin 
zone as if hole pockets exist. Each (fictitious) Fermi surface is discretized into 16 patches as in 
ref. |l6llll[l4l. In the calculation of RG flow these hole band patches represented by the fictitious 
Fermi surface points are treated in the same way as electron band patches. In the second stage 
where A < AEq the hole bands no longer contribute virtual states in the one-loop diagrams. 
When we compute pairing the fictitious Fermi surface are dropped. Thus in the second stage of 
the FRG the flow of pairing channel is essentially the marginally relevant Cooper channel flow 
discussed by Shankar [|3T1l . 

In our previous FRG study for systems involving both electron and hole pockets, the final 
effective interaction is captured reasonably well by a Ji — J2 model [fT4l| (as far as the antiferro- 
magnetism and pairing are concerned). In addition, in those cases, the final effective interaction 
is a result of the FRG flow over the entire bandwidth, not just in slivers around the Fermi 
surfaces. Thus to gain a qualitative understanding of the pairing in K^.Fe2-j;Se2 one can use 
an effective Ji — J2 interaction fTl to replace the FRG flow for cutoff A > AEq. We have 
analyzed the pairing form factor using such a model and find that for Ji, J2 that are both anti- 
ferromagnetic and weak, the nodeless c/^.2_^2 is the dominant pairing channel for the Ka;Fe2Se2 
bandstructure as long as J1/J2 > 0.213. In contrast the s± is dominant in the LaFeAsO band 
structure from J1/J2 = up to J1/J2 ~ 1. Of course, in general what affects the final pairing 
strength and form factor involves both the strength and form of the effective interaction as well 
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Figure 3: (Color online) FRG flow of two leading pairing channels for (f/i, U2, Jh)={^, 2, 0.7) 
eV. Horizontal axis is the negative of the decadic logarithm of RG cutoff A. Vertical axis is the 
effective pairing strength (c.f. ref. [il4ll ). 

as the geometry of the Fermi surfaces. 

4 FRG Results 

The presentation of our results follows the format in our previous publications |l6l HI [IS. For this 
FRG calculation with the parameters and RG scheme described above, we obtained divergent 
flow of several pairing channels, the flow of the two leading channels are plotted in fig.[3l The 
strongest pairing channel has the d^^-yZ symmetry. The associated divergence energy scale is 
~ 4meV. The second strongest pairing is a extended s-wave with form factor approximately 
described by cos{kx) + cos{ky). It has four nodes on each electron Fermi surface. The pairing 
form factors of both channels are depicted in fig. IH Interestingly both pairing channels are 
favored by the nearest-neighbor antiferromagnetic interactions Ji. In fact as long as J1/J2 > 
0.2 the Ji — J2 effective interaction predicts a leading dx2_y2 pairing whose form factor is almost 
identical to our FRG result. The overlap between the form factor predicted by the Ji — J2 model 
and that determined by FRG is shown in fig. [51 
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Figure 4: (Color online) Gap form factors (c.f. ref. (14]) of the two leading pairing channel 
for {Ui, U2, o/h)=(4, 2, 0.7) eV. Horizontal axis 9 is the polar angle on each Fermi surface with 
respect to its center. 9 = indicates the +kx direction. Upward(red) and downward(green) tri- 
angles label the two electron Fermi surfaces at (0, tt) and (tt, 0) respectively. Top: The Leading 
dx2_y2 gap. Bottom: The subleading extended s-wave gap. 




Figure 5: (Color online) The overlap between the gap form factors of leading pairing channels 
obtained by FRG [{Ui, U2, Jh)=(4, 2, 0.7) eV] and Ji — J2 mean field theory versus the ratio 
J1/J2 (J2 > is antiferromagnetic). 
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5 Discussions 



There are evidences from NMR that the antiferromagnetic fluctuation is weak in KJe9_,^Se9 [l32ll 
On the surface this seems to be at odds with the notion that spin fluctuation is the main source 
of pairing. However, NMR probes spin fluctuations at very low energy. It is the high energy 
spin fluctuations that are good for pairing. 

Obvious implications of the nodeless dx2^y2 pairing is the presence of (tt, tt) neutron res- 
onance in the superconducting state. In addition, phase sensitive measurement utilizing the 
relative orientation of crystallographic grains, or the detection of randomly trapped half flux 
quanta in polycrystalline materials, can be done similar to the cuprates. Fourier transform scan- 
ning tunneling spectroscopy should in principle reveal the enhancement of the (tt, tt) scattering 
by the magnetic field. Whether any of these is true is remained to be studied. 

Finally we would like to list a few caveats. First, there is no a pnor? justification for ignoring 
the (0, 0, tt) electron pocket found in the DFT calculation. It is important to study the effect of 
the three dimensional electron pocket on the pairing symmetry. Second, although the particular 
parameter set (including the chemical potential) we used in this study yielded drc2_y2 pairing, 
it is important to investigate the robustness of this pairing symmetry against chemical potential 
shift and/or the changing of interaction parameters. Presently we have studied a different set of 
interaction parameters {Ui,U2, Jh)=(4, 4, 0) eV and found a weaker pairing divergence (fig.©. 
More importantly the leading pairing symmetry becomes a fully gapped s-wave, and the sub- 
leading one becomes dxy (fig- 12]). Both these channels may be driven by strong J2. The fully 
gapped s-wave can be viewed as the remnant of the s± pairing after the central hole pockets 
are removed. This result clearly calls for a more systematic study of the dependence of pairing 
symmetry on the parameters in our model. Third, the real materials have Fe vacancies, and even 
possible partial Fe vacancy ordering. This is obviously ignored in the current study. Last, the 
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Figure 6: (Color online) FRG flow of two leading pairing channels for {Ui, U2, Jh)=(4, 4, 0) 
eV. Horizontal axis is the negative of the decadic logarithm of RG cutoff A. Vertical axis is the 
effective pairing strength (c.f. ref. [il4in . 

effects of impurity scattering on destabilizing the d^i-y^ pairing clearly need to be addressed. 

In summary we have studied the pairing instabilities in K2.Fe2-ySe2 using a two stage func- 
tional renormalization group method. Our calculation suggests the leading pairing channel is 
dx2-y2 and the subleading one is extended s. We find the buried hole bands make important 
contributions to the final effective interaction because of the closeness of their maxima to the 
Fermi energy. The leading d^i^yi pairing form factor is well captured by an effective Ji — J2 
interaction with dominant nearest neighbor antiferromagnetic exchange. At the end of the paper 
we briefly mentioned several obvious experiments that can test the dx2_y2 symmetry. We also 
list some caveats of the present study. 

DHL is supported by DOE grant number DE-AC02-05CH1 1231. FY is supported by the 
NSFC Grant No. 10704008. ZYL and TX are supported by National Natural Science Foundation 
of China and by National Program for Basic Research of MOST, China. 
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